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ABSTRACT 
Surface  wave  action  has  been  proven  to  be  the  biggest  deterrent 
to  a  ship's  progress  through  the  water.   Until  the  present  time  fore- 
casts of  surface  waves  have  been  available  for  only  a  few  days.   In 
order  to  determine  an  optimum  ship  route  across  the  ocean  these  wave 
forecasts  must  be  extended  to  15-20  days.   Using  three  computer  programs 
from  the  Fleet  Numerical  Weather  Facility  and  the  U.S.  Weather  Bureau 
5-day  and  30-day  forecasts  of  the  mean  sea  level  pressures  as  input 
data,  surface  wave  forecasts  may  be  prepared  for  periods  up  to  30  days. 
When  compared  with  the  actual  waves  over  the  specified  period  of  time, 
the  forecast  wave  patterns  show  more  skill  than  forecast  waves  derived 
strictly  from  climatology. 
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1.   INTRODUCTION 

Optimum  track   ship  routing    (OTSR)    has   been   in  exist ance   in  the 
U.S.    Navy   for  almost   ten  years    under  this    specific   name.      An  optimum 
track   is  one  chosen  from  port   to   port   across   the  ocean  that   will   allow 
a  ship  to  traverse   the  route   in  the   least   amount  of  time  and   also   en- 
counter  less   unfavorable  weather.      James  [_5j  has   shown  that   a  definite 
increase   in  operational   efficiency   is  made  possible   by  the  routing  of 
transoceanic   vessels   along  tracks   based  primarily  on   long-range  ocean 
wave  forecasts   rather  than  along   standard   seasonal   tracks.      Such  predic- 
ted tracks   are  believed  to   be   superior  to  the   standard   seasonal   routes 
in  terms  of  travel   time,    passenger  comfort,    damage   to   cargo   and   ship, 
and   fuel   consumption.      Absolute  values   for  the  reduction  of  the  above 
items    are  difficult   to   present,   but   certainly  the  great   weight  of  evi- 
dence presented   in  his   evaluation  supports   the  contention  that   efficiency 
of  ship  operation   is    improved  by  oceanographic   ship  routing.      The   Mili- 
tary Sea  Transportation    Service  tentatively  estimates   that   each  hour 
saved    in  ship  operations   is   a   savings  of   $200.00  to   them. 

Some  of  the   international    steamer  tracks    in  use  today  are  modifi- 
cations  of  routes    laid  out  by  Lt .    Matthew  Fontaine  Maury,   USN,    as   far 
back  as   the   1850's.      Maury's   routes   were  based  on  ships'    logs,   clima- 
tology and  statistics.      Today,    James  \_^j  has    found  that   the  most    impor- 
tant  parameter  retarding  a  ship's   progress   is   surface  wave  action. 

Wave   action  actually  affects   a   ship's   progress    in  two   ways.      First 
simply  by  retarding  the  ship's    speed   through  the  disturbed  water  and 
secondly  by  making  the  Captain  order  a  reduction   in  speed  due  to  the 
violent  effects  of   roll,    pitch,    and  other  associated  motion  of  the   ship 
in  creating   eddy  resistance. 


Ships   not  utilizing  OTSR,   when  making  a  transoceanic  voyage  in  the 
Northern  Hemisphere,   usually  travel  a  more  northerly  route  in  the  summer 
months   and  then  shift  to   a  southerly  route  in  the  stormy  winter  months. 
This  usually  results   in  avoiding  the  largest   storms  but   at  a  great  cost 
in  time  and  money.     Reduction  in  weather  damage  is  not  the  whole  solu- 
tion, however,   since  the  ideal   situation  is  to  save  time  and  increase 
the  distance  traveled  per  day.     Often  times,   the  weather  situation  may 
be  such  that  the  off-season  route  has  the  best  sea  conditions. 

At  the  present  time  a  5-day  forecast  of  the  surface  weather  phenom- 
enon is  the  maximum  utilized  in  ship  routing.     If  this  could  be  extended 
say  to   15-20  days,   an  entire  voyage  could  be  routed  across  the  Pacific 
Ocean  to  avoid  storm  areas  and  reduce  transient  time.     Operationally, 
the  procedures  for  ship  routing  are  slightly  different  in  the  Atlantic 
and  Pacific  operating  areas.     More  emphasis  is  placed  on  storm  track 
trends  and  climatic  means  in  the  Pacific  since  the  5-day  forecast  covers 
less  than  half  the  normal  crossing. 

It  is  the  purpose  of  this  paper  to  derive  extended  range  forecasts 
of  wave  heights,  directions  and  periods  for  15-20  days  in  advance  of  a 
proposed  sailing  across  the  Pacific  Ocean. 

Naturally,  constant  surveillance  would  have  to  be  maintained  over 
the  routed  vessels  as  the  5-day  forecasts  are  up-dated  and  when  sea 
conditions  warrant,   the  routes  would  be  modified  to  indicate  a  superior 
route. 


2.      COLLECTION  OF  DATA 

The   U.S.    Weather  Bureau  publishes   a   5-day   forecast  of  the  mean   sea 
level   pressures   around  the  Northern  Hemisphere  and  this    is   presented  on 
individual   maps   for  each  of  the   five  days.      These  maps   are  distributed 
on  Monday,    Wednesday,    and  Friday  of   each  week   and   are   portrayed  on  a 
polar  stereographic   1:50,000,000   scale.      The   forecasts   are  based  on  a 
mean  circulation  method  by  averaging  ten  twice-daily  values   of   pressure 
at   standard   intersections  of   latitude   and   longitude  on   synoptic   sea 
level   weather  maps.   |_1  \  The  Weather  Bureau  also  makes    30-day  forecasts 
using  methods   similar  to   those  employed   in  preparing  the   5-day   fore- 
casts. 

These  maps   were  utilized  to  obtain  a   surface  wave  forecast   for  an 
extended   period  of   time   across   the  Pacific  Ocean.      The  period  of   time 
selected  was   26  July   1966  to   06  August    1966,    a  total   of    12   days.      The 
U.S.    Navy  Fleet   Numerical  Weather  Facility   located   at   Monterey,   Califor- 
nia  publishes   24  and  48-hour  wind-wave   prognoses   which  were   used   for  the 
first   two   days  of  the  forecast,    26   and  27   July   1966.   {_2^  For  28,    29,    and 
30   July   1966  the   last   three  maps   of   the  U.S.    Weather  Bureau  5-day   fore- 
cast  series  were  applied.      For  the   31st  of  July  through   6  August   the 
30-day  mean   sea   level    prognostic   chart  dated    15  July   1966  was   utilized 
on  a   day   to   day   basis. 

In  the  following  explanation,    focus    is   placed  only  on  one  map   in 
showing  how  the  surface  wave   forecast  was   derived.      The   remaining  days 
of   the  forecast   were  treated   in  a  similar  manner. 

The   first   step   in  arriving   at   the   forecast   was   to   use  the  FNWF 
FAREB  program.      This   program  was   written  by  Harry  N.    Farnsworth  on   1 
May   1965   at   the   Fleet   Numerical   Weather  Facility.      The   basic   purpose 


of  this  program  is  to  reconstruct  any  scalar  field  in  memory  in  a  CDC 
1604  computer  as  a  3969-point  field.  An  attempt  will  not  be  made  to 
explain  each  step  of  this  program  in  detail,  however,  input  and  output 
values  along  with  an  example  for  preparation  of  a  field  will  be  dis- 
cussed. 

Input  values  were  extracted  from  the  source  chart  by  employing  a 
plastic  overlay  of  the  standard  FNWF  63  x  63  grid.   If  a  grid  overlay 
were  not  available  latitude-longitude  input  could  be  used.  A  proper 
interval  should  be  selected  in  order  to  combine  sufficient  pattern  de- 
lineation with  a  minimum  of  data  inputs.  The  initial  guess  field  is 
considered  in  the  interval  selection  and  since  a  guess  field  of  zero 
was  used,  an  interval  of  two  grid  points  or  every  other  point  was  chosen. 
Then  by  superimposing  the  grid  overlay  on  the  source  chart  an  initial 
grid  point  was  established  and  grid  point  values  were  read  along  a  J-row 
at  a  specified  point  spacing  (every  other  point).  A  new  initial  point 
may  be  established  whenever  desired  to  reorient  the  input  data.  The 
input  points  selected  were  denoted  by  dots  and  dashes  on  the  FNWF  over- 
lay.  Since  the  source  chart  was  in  millibars,  the  values  were  read  to 
the  nearest  millibar  interpolating  where  necessary.  When  all  values  had 
been  read  from  the  source  chart  in  the  desired  area  and  tabulated  they 
were  converted  to  1000-mb  heights  using  a  standard  conversion  of  7.5  mb 
per  200  feet  with  a  reference  point  of  1000  mb  equal  to  zero  feet.  All 
values  had  to  be  coded  as  four  digit  numbers  for  entry  on  the  standard 
IBM  card.  The  values  were  assumed  to  be  positive  unless  preceded  by  a 
minus  sign.  There  were  some  low  centers  present  with  central  pressures 
below  1000  mb  and  these  were  coded  accordingly. 

In  order  to  assemble  the  card  deck  the  identification  and  scaling 


information  was  prepared.   On  the  initial  card  the  titling  information 
and  initial  guess  value  were  placed.   The  standard  procedure  is  for  the 
titling  information  to  consist  of  the  date-time  group,  some  identifying 
number,  in  this  case  D1000,  and  then  any  explanatory  information  that 
might  be  desired  such  as  name,  etc. 

The  program  was  then  put  into  operation  and  the  final  field  con- 
structed through  smoothing,  relaxation,  and  a  combination  of  the  two. 
A  grid  print  in  feet  of  the  height  of  the  1000-mb  surface  was  produced 
and  the  field  packed  and  written  on  magnetic  tape.   As  mentioned  pre- 
viously the  main  objective  of  this  program  was  to  get  the  data  on  tape 
in  an  organized  field  in  order  to  go  on  to  the  next  step. 


3.   WIND  EXTRACTION 

The  next  step  in  arriving  at  the  wave  forecast  was  to  utilize  the 
FNWF  wind  extraction  program.   By  using  the  1000-mb  constant  pressure 
surface  as  the  input  this  program  computes  the  geostrophic  u,  v  wind  com- 
ponents, as  well  as  the  direction  and  speed  of  the  wind  by  solution  of 
the  geostrophic  wind  equation.  There  are,  however,  some  limitations  to 
this  program.  Wind  speeds  that  are  computed  by  the  geostrophic  wind 
equation  are  generally  greater  than  actual  speeds  in  regions  of  cyclonic 
trajectory  curvature  and  less  than  actual  speeds  in  regions  of  anticy- 
clonic  curvature. 

This  program  was  run  on  the  3200  computer  at  the  Fleet  Numerical 
Weather  Facility.  An  example  of  a  3200  run  sheet  is  shown  in  Appendix 
I.   Some  modifications  had  to  be  made  in  order  to  read  the  mean  D1000 
field  and  the  climato logical  D1000  field  taken  from  the  U.S.  Weather 
Bureau  records.   These  modifications  were  made  on  a  3200  octal  patch 
and  an  example  of  this  is  shown  in  Appendix  II.  The  output  of  this 
program  was  packed  and  written  on  magnetic  tape  for  an  input  to  the  wind 
waves  program. 
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4.   WIND  WAVES 

The  next  step  was  to  derive  the  final  product.   The  Fleet  Numerical 
Weather  Facility  Wind  Waves  program  with  slight  modifications  was  run  on 
the  Control  Data  Corporation  1604  computer  to  accomplish  this  purpose. 
In  order  to  eliminate  some  of  the  steps  not  needed  in  this  program  a 
flex-tape  patch  (Appendix  III)  was  made  and  inserted  into  the  computer 
at  the  appropriate  time. 

The  program  input  consisted  of  three  magnetic  tapes.   Tape  number 
one  contained  the  1000-mb  field  obtained  in  the  FAREB  program  previously- 
mentioned.   Tape  number  two  contained  the  climatological  sea  surface 
temperatures  for  the  month  of  July.   Tape  number  three  was  comprised  of 
the  u  and  v  components  of  the  wind.   A  fourth  tape,  blank  pool  tape, 
was  used  in  order  to  copy  the  output  of  the  program.   This  was  necessary 
in  making  a  difference  field  with  the  actual  analyzed  waves  in  order  to 
determine  the  merits  of  the  5  and  30-day  forecasts  for  preparing  long 
range  wave  prognoses. 

The  Waves  program  corrects  the  geostrophic  wind  directions  for 
surface  friction  by  assuming  a  constant  cross-contour  flow  of  15  de- 
grees toward  lower  pressure.   A  stability  correction  is  then  made  by 
reducing  the  wind  speed  depending  upon  the  general  air  mass  type  and 
the  sign  of  the  surface  temperature  advection.   For  example,  a  mass  of 
maritime  tropical  air  blowing  toward  colder  sea  surface  temperatures 
would  warrant  a  greater  reduction  in  wind  speed  than  a  mass  of  maritime 
polar  air  blowing  toward  warmer  sea  surface  temperatures.   The  duration 
of  the  wind  is  determined  at  each  grid  point  by  computing  the  length  of 
time  the  wind  has  been  blowing  from  the  same  direction  (within  t   22  de- 
grees).  The  maximum  duration  has  been  set  at  18  hours  according  to  the 
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French  school  of  thought  which  considers  all  seas  to  be  fully  developed 
after  this  length  of  time.  An  effective  wind  speed  is  then  obtained  by 
an  averaging  technique.  The  significant  wave  height  and  wave  period  is 
computed  using  two  singular  relationships  established  at  the  Fleet  Nu- 
merical Weather  Facility: 

Hl/3  =  K1v2d  +  K2V  (4,1) 

Tl/3  =  V(K3  +  K4D)  +  K5  (4,2) 

where  HW3  and  T\  ,*   are  tne  significant  wave  height  and  period,  respec- 
tively; V  is  the  effective  geostrophic  wind  speed;  D  is  the  duration 
time  and  Kn  to  Ke  are  empirical  constants  tuned  to  reduce  the  root  mean 
square  error. 

An  example  of  a  1604  run  sheet  is  shown  in  Appendix  IV.  The  tape 
set-up  and  instructions  listed  were  actually  used  to  run  this  program. 
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5 .   VERIFICATION 

It  was  decided  at  the  onset  that  once  the  forecast  wave  pattern  was 
derived  a  comparison  would  be  made  of  the  derived  field  and  the  actual 
analyzed  waves  for  the  specified  period  of  time.   To  accomplish  this 
purpose  an  A-B  program,  called  SAK  and  written  by  R.  McDonald  at  Fleet 
Numerical  in  1962,  was  utilized.   This  program  computes  a  difference 
field  using  any  two  given  fields,  A  and  B,  a  pillow  and  root  mean  square 
error  for  the  entire  grid.   Pillow  and  RMSE  are  defined  as  follows: 

x 

Pillow  =  ^_      (A  -  B)n         x  =  Dependent  on  number   (5.1) 

".  of  points  differenced 

n=l  r 


£      [(A-B)  -  Pillow!  (5.2) 


n=l 


RMSE  = 

\ 

A  further  comparison  was  made  of  the  actual  analyzed  waves  and  a 
wave  forecast  based  solely  on  climatology.   Sea  level  pressure  normals 
were  obtained  from  the  U.S.  Weather  Bureau  and  a  D1000  field  was  de- 
rived.  From  this  field  a  wave  forecast  was  made  following  the  same 
basic  steps  as  explained  above. 
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Pacific  Ocean  region  and  the  rest  of  the  grid  left  blank.  When  the  SAK 
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APPENDIX  I 

3200  Run  Sheet 
Run  #     12      Date   9-26-66    Est.  Time  10   Stop  if  exceeds 

Programmer Truax Actual  Time 

Rush Routine Rockpile Listing  Only  

(  )  BPI IN OUT List  &  Binary  Deck  

(  )  BPI IN OUT List  &  Mini  Tape  


(  )  BPI    IN    OUT        Stack  for  Transmit 


(    )  BPI  IN OUT 


Program  Requires:   Drum       Printer 


1.     Load:      Binary  Tape Binary  Deck  Patch  Cards 


2.     Load:    Drum  Records 

Tape  Lib  Records 


3.      I&C  Loc =  ;   Loc_ 

Loc -  ;  Loc 


4.      Set  JK  :    SK      :    Other 


5.      Set  Wds    13,l4(Dtg) Wdl5 Wdl6 Wdl7_ 

Set  A  =  Q  =  _____ 


Bl   =  B2  =  B3  = 


BPI   =  BPO  =  STO  = 


Start  P  ■  Expected  Final Dump  Yes  No 

Dump  to Mode ;  to Mode 


Other  Instructions: 

(1)  Set  up  for   "FWBXB"  using  tape  H-24  (200  Density)   as   input, 

(2)  Read  in  Bi-Octal  correction  cards. 

(3)  Run  wind  31A  on  output  fields  from  FWBXB. 

(4)  Save  output  on  enclosed  ST-500  tape. 
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APPENDIX  II 
3200  OCTAL  PATCH  FOR  WIND  EXTRACTION 


02435 
02436 
02437 
02440 
0  2441 
02442 
02443 
02444 
02445 
02446 
02447 
02450 
02451 
02452 
02453 
02454 
01007 


64446561 
45202020 
00000170 
00000036 
00567000 
00160400 
63121220 
11005305 
64634371 
44612320 
00000170 
00000036 
00567000 
00160400 
63121220 
11005305 
00010201 
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APPENDIX  III 


1604  FLEX  PATCH  FOR  WIND  WAVES 
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APPENDIX  IV 

1604  Run  Sheet 
Run   #7  Date  9-16-66   Est.  Time  10    Stop  if  exceeds 

Programmer     Truax Actual  Time  

Rush Rout  ine   X   Rockpile Assembly:  AR Scrap 


Ch  3/4  Bioct. Bin. Cards 

Ring 
(1)   1000  mb   BPI  200  in.      out  X _ 

Ch  5/6 

Ring 

( 3)  Pool  BPI  556  in.  X  out 

(4)  U,V  BPI  200  in. out  X 


1.  Load  Drum  library  records  "WAVES"  (Operational  program) 

2.  Load  bioctal  tape(s) Flex  tape(s)   „one„, 

3.  I&C  Location       A(L)  A(R) 


(2) 

SST 

BPI  200 
BPI 
BPI 

in. 
in. 
in. 

out 
out 
out 

X 

(  ) 

(  ) 

A(L) A(R) 


4.  Set  Wd.  10  =  Wd.  16  Pos Neg_ 

5.  Set  Breakpoint JK  1     Jl 3__  SK  __  _ 

6.  Set   A  Set  Q 


7.  Start  P  *  1026  Expected  Final  616   Dump  Yes  X   No_ 

8.  Dump   600    to    10000    Instr/Data 


Other  Instructions: 

1,  Start  at  1026  program  will  stop  at  616,  hit  go 

2.  Save  pool  tape  if  program  runs 
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APPENDIX  V 
COMPARISON  OF  RMSE  &  PILLOW 


Date 


Derived  Waves 


Climatological  Waves 


RMSE  (FT)    PILLOW  (FT)     RMSE  (FT)    PILLOW  (FT) 


28  July 

29  July 

30  July 

31  July 

1  August 

2  August 

3  August 

4  August 

5  August 

6  August 


3 
3 
3 
3 
3 
3 
3 
3 
3 
3 
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